Turbulence-induced magnetic flux asymmetry at nanoscale junctions 
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It was recently predicted [J. Phys.: Condens. Matter 18, 11059 (2006)] that turbulence of the 
electron flow may develop at nonadiabatic nanoscale junctions under appropriate conditions. Here 
we show that such an effect leads to an asymmetric current-induced magnetic field on the two sides 
of an otherwise symmetric junction. We propose that by measuring the fluxes ensuing from these 
fields across two surfaces placed at the two sides of the junction would provide direct and noninvasive 
evidence of the transition from laminar to turbulent electron flow. The flux asymmetry is predicted 
to flrst increase, reach a maximum and then decrease with increasing current, i.e. with increasing 
amount of turbulence. 

PACS numbers: 73.23.-b, 47.27.Cn, 73.63.Rt 
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The hydrodynamics of the electron liquid dates back 
to earlier studies by Madelung, Bloch P, Q and later 
on by Martin and Schwinger [3|. In this latter work, 
in particular, it was shown that the many-body time- 
dependent Schrodinger equation (TDSE) can be written 
exactly in hydrodynamic form in terms of the density 
n{r,t) and velocity field v(r,t) = j(r, i)/n(r, i), the ratio 
of the current and charge density, with all many-body 
interactions lumped into a two-particle stress tensor. 

In recent years, the analogy of the electron flow 
with classical fluid dynamics has been pushed even fur- 
ther with the development of time-dependent density- 
functional methods and the consequent realization that 
under certain conditions, the exchange-correlation poten- 
tials can be written in hydrodynamic form [3, fsj . More 
recently, it was shown that electron flow in nanoscale con- 
strictions satisfies the conditions to write the two-particle 
stress tensor in a form similar to the stress tensor of the 
Navier-Stokes equations with an effective viscosity of the 
electron liquid (see also below) 0. The most striking 
prediction of this result is that, under specific conditions 
on the current, density and junction geometry, the elec- 
tron flow should undergo a transition from laminar to 
turbulent regimes, [ij Recently, this behavior was con- 
firmed numerically by solving directly the TDSE within 
time-dependent current-density functional theory [tI and 
comparing the results with the generalized Navier-Stokes 
equations derived in Rcf. 6]. In experiments, however, 
detecting turbulence via direct imaging of the current 
density remains challenging. For instance, scanning- 
probe microscopy (SPM) experiments which image the 
current flow in a 2D electron gas (2DEG) have been re- 
ported [8[. These experiments employ an SPM tip to 
reflect electrons back toward the junction, and measure 
the resultant change in the total current. This means 
that the image thus obtained gives the correlation be- 
tween the tip position and junction current, which docs 
not necessarily correspond to the magnitude of the cur- 
rent density. Moreover, SPM-type experiments necessar- 



ily disturb the electron flow at the tip position, and are 
therefore essentially invasive. 

Another way to probe turbulence would be to measure 
the noise properties of the current. However, the analyt- 
ical description of turbulence is a notoriously intractable 
problem, thus making it unclear what noise character- 
istics the turbulent flow of electrons would generate. In 
addition, this noise would necessarily correlate with other 
intrinsic types of noise, especially shot noise 

In the present Letter we show that the measurement of 
the current-induced magnetic field at the two sides of an 
otherwise symmetric nano junction provides a direct and 
non-invasive way of measuring the transition from lami- 
nar to turbulent flow. In particular, we predict that the 
fluxes ensuing from the current-induced magnetic field 
across two surfaces on the two sides of the junction would 
at first become increasingly different with increasing cur- 
rent. This asymmetry reaches a maximum, and then 
decreases with further increase of the current. The mea- 
surement of these fluxes is within reach of present exper- 
imental capabilities, and thus the observation and study 
of this phenomenon would provide valuable insight into 
the transport properties of nanoscale systems. 

The structure we have in mind consists of two symmet- 
ric regions of a 2DEG connected non-adiabatically by a 
nanojunction (the edges of this structure are represented 
with solid lines in the left panel of Figs. [T]and[2|). The 
non-adiabaticity requirement is due to the fact that, as 
shown in Ref. [6|, an adiabatic constriction produces a 
Poiseuille flow, which is laminar for essentially all cur- 
rents one can effectively inject in a 2DEG. 22] The lat- 
eral (y-direction) boundaries are closed to current flow, 
and the longitudinal (cc-direction) boundaries are open, 
with current being injected in the "top" boundary and 
exiting in the "bottom" boundary. We then envision two 
identical surfaces - placed at a given distance from the 
2DEG in the z direction - across which we calculate the 
current-induced magnetic flux (see Fig. [1]). These mag- 
netic fluxes can be measured by two superconducting in- 
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terference devices (SQUIDs) fio'l located on the two sides 
of the junction as illustrated in Fig. [TJ 

Our starting point is the time-dependent Schrodinger 
equation written in the approximate Navier-Stokes form 
for an incompressible fluid |6i] 

An(r,t) =0, V-v(r,t) = 0, 
d 

m*n(r, t)DtVi{Y, t) ^ - — P(r, t) + t]V'^v,{y, t) 
- n{r, t) — Vcxt{r,t) 

where Dt = ^ + (v • V) is the convective derivative, 
m* is the effective mass of the electrons, n(r, t) is the 
electron density, P(r, t) is the pressure of the electron 
liquid, and Kxt (r, t) is the sum of the Hartree and the 
ionic potentials. The quantity r] = Hnf{n) is the viscos- 
ity of the electron liquid, with f{n) a smooth function of 
the density. The values of the viscosity as a function of 
density have been calculated using linear-response the- 
ory llll.[T3 |: here, we use the 2D interpolation formula of 
Ref. We also employ the jellium approximation for 
the electron liquid, which together with the assumption 
of incompressibility, allows us to neglect spatial varia- 
tions of T4xt(r,i)- Incompressibility of the electron liq- 
uid represents to a good approximation the behavior of 
metallic quantum point contacts 0| [l^l • 

The current density was calculated numerically as a 
solution of Eqs. H]) |13| for a nanojunction 28 nm wide. 
We have used Dirichlct boundary conditions for the ve- 
locity at the inlet, and Neumann boundary conditions at 
the outlet. We use parameters corresponding to a GaAs- 
based 2DEG: m* = 0.067me, and n = 5.13 x 10" cm'^. 
The calculations were performed at a fixed value of the 
total average current flowing through the system selected 
in the range from 0.001 iiA to 10 /LtA. The magnetic field 
profile was found for each calculated current density dis- 
tribution. The size of the surface area across which we 
calculated the magnetic flux was chosen to be 200 x 200 
(nm)^. Each surface is displaced laterally to one side of 
the nanojunction as shown in Fig. [T] This surface rep- 
resents the SQ UID area and is within reach of current 
technology [14|, although, in principle, the SQUID's area 
does not necessarily need to be so small. We assume that 
the surfaces are located 50 nm above the 2DEG; there- 
fore, the distributions of magnetic field were calculated 
at this distance from the 2DEG. Increasing this distance 
by a factor of two decreases the flux by about a factor of 
two. [2J] For convenience, from now on we shall call these 
two surfaces, SQUID 1 area and SQUID 2 area (Figs. [1] 
and [2]). 

The magnetic field fluxes allow us to characterize the 
degree of asymmetry in the current flow pattern as well as 
to probe such specific features of turbulent current fiow 
as eddies. We find that at low currents, the electron fiow 
pattern is symmetric (see left panel of Fig [T]) . More pre- 
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FIG. 1: (Color online) Laminar current flow at a low value 
of the total current (0.1 /xA). Left panel: electron velocity 
distribution. The arrow length is proportional to the velocity 
magnitude. Right panel: normal component of the magnetic 
field through a plane 50 nm above the 2DEG. The surface 
areas indicated with SQUID 1 and SQUID 2 represent two 
areas across which we calculate the magnetic flux (see text). 

cisely, the patterns of current flow in the two electrodes 
are mirror images of each other, with the overall sign of 
the fiow reversed This is illustrated in the left panel 
of Fig. [U where we plot the velocity distribution of the 
electron liquid at a simulation time of 24 ps. The total 
current is 0.1 /iA, which is small enough that the flow is 
in the laminar regime, as it is evident from the figure. In 
the right panel of Fig. [1] we show the z-component of the 
magnetic field through a plane 50 nm above the 2DEG. 
As is typical of the laminar regime, the current-induced 
magnetic fields above the top and bottom contacts are 
almost symmetric with respect to the center of the junc- 
tion producing an almost symmetric flux across the areas. 
Hence, SQUIDs positioned as indicated in Fig. [1] would 
measure almost equal magnetic fluxes. 

By increasing the current, the current density in the 
source and drain sides loses top-down symmetry: the 
current density in the outgoing side becomes turbulent, 
while the current density in the incident side remains 
laminar (Fig ^ . We note that at large currents we ob- 
serve the formation of turbulent "eddies" which evolve 
in time, rather than a completely chaotic current density 
distribution. This means that at the current values we 
consider here turbulence is not fully developed [t!]. 

Fig. [2] illustrates the behavior typical for the turbulent 
regime. This plot corresponds to a total current of 1.0 fiA 
at time t = 2A ps. As before, the left panel illustrates the 
electron velocity distribution. Unlike the electron veloc- 
ity field presented in Fig. [U the electron velocity distri- 
butions in the top and bottom electrodes in Fig. [2] are 
no longer symmetric. In particular, the electron velocity 
distribution in the bottom contact shows eddies and an 
increased current density in the middle of the junction. 
Such a velocity field is responsible for a much stronger 
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FIG. 2: (Color online) Turbulent current flow at high value of 
the total current (l^A). Left panel: electron velocity distri- 
bution. Right panel: normal component of the magnetic field 
through a plane 50 nm above 2DEG. The top-down symme- 
try in magnetic field distributions across the SQUID areas is 
lost. The magnetic field flux through SQUID 2 is significantly 
higher then the flux through SQUID 1. 



magnetic field above the bottom contact (in particular, 
in the SQUID 2 area in Fig. The magnetic field dis- 
tribution above the top contact in Fig.[2]is "smooth" and 
uniform, and has a structure similar to the magnetic field 
distribution above the top contact of Fig. [TJ 

For both the laminar and turbulent cases, as time 
passes, the fluxes through the top and bottom SQUID 
areas saturate to constant values. This is shown in the 



5 I ..I I 1 I 




Time (ps) 




+T-r| 1 I 1 I 1 I ' I 

10"' 10"' 10"' 10° lO' 

/(^A) 

FIG. 3: (Color online) Asymptotic value of the ratio of the 
flux through SQUID 2 area to the flux through the SQUID 1 
area as a function of the total current. See text for details. 
Inset: Time dependence of fluxes (in units of <l>o = h/2e) 
through the top (dashed blue line) and the bottom (solid red 
line) SQUID areas, for a total current of 1.0 /iA. 
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FIG. 4: (Color online) Magnetic fleld correlation function, 
for total currents of 1.0 fi A (solid orange curve) and 0.1 /lA 
(dashed blue curve). Because the total current differs between 
the two by a factor of 10, we have scaled the 0.1 ^A curve by 
a factor of 10^. Even after this scaling, Czz is signiflcantly 
larger in the turbulent (high-current) case than it is in the 
laminar (low-current) case. 



inset of Fig. [3] where we plot the flux through the top 
and bottom SQUID areas as a function of time for the 
1.0 A case. We can now determine the magnetic signa- 
ture of the transition between the laminar and turbulent 
regimes by plotting the asymptotic value of the ratio of 
the magnetic fluxes through the two SQUID areas as a 
function of the total current. This is illustrated in Fig. [3] 
and shows the main flndings of this paper. 

At low currents (laminar regime), the current den- 
sity in the top and bottom contacts is highly syinmetric. 
Therefore, the ratio of the fluxes is near unity (25|. There 
is a definite transition at the critical current Ic i:: 0.3 /iA, 
after which the ratio of the magnetic fluxes increases by 
more than 300%. The flux ratio reaches a maximum 
at a current of about IfiA. Increasing the current fur- 
ther, we flnd that the flux ratio decreases. This can 
be explained as follows. Near the critical current, the 
eddies are small and stationary. These smaller, local- 
ized eddies lead to a higher ratio <i>2/$i and the in- 
crease we observe. On the other hand, by increasing 
the current further, the eddies grow in size and "spread 
out" away from the junction thus reducing the flux ra- 
tio. Note that fully-developed turbulence would decrease 
the ratio ^2/^1 even further, since eddies having clock- 
wise and counter-clockwise senses would be continuously 
generated and destroyed. In addition, the noise level of 
modern SQUIDs [ll is below lO^^^o ( = h/2e is 
the magnetic flux quantum), which is well below the flux 
magnitude shown in the inset of Fig. [3l 

We conclude by quantifying the degree of turbulence of 
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the current-induced magnetic field. While this cannot be 
directly measured, it provides insight into the properties 
of the turbulent regime attainable experimentally. Let 
us then calculate the magnetic field correlation tensor, 
which quantifies the spatial correlation of the magnetic 
field at different points in space. We define this tensor as 
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{{B,ir) - B,{r + <5r))(S,(r) - B,{r + Sr))) . (2) 



Here, Bi and Bj denote components of the magnetic field, 
and 6r is a given vector. The brackets (. . .) denote aver- 
aging over all pairs of positions separated by Sr within a 
given region. Note that, even before performing a spatial 
average, the magnetic field already has a nonlocal char- 
acter, in that the magnetic field at a point is due to the 
velocity of charges in the whole system. 

In Fig. m we plot the magnetic field correlation ten- 
sor Czz, for i5r = (fe,0) at 50 nm above the 2DEG as 
a function of Sx. The spatial averaging was carried out 
in a region "downstream" from the junction, in the re- 
gion X = [— 487nm, — 235 nm], y = [0nm,259nm]. As 
expected, in the laminar case, the magnetic field varies 
with distance by a small amount, so that Czz is small. 
In the turbulent case, the presence of the eddies leads to 
a magnetic field that correlates spatially, causing Czz to 
increase with distance. 

Finally, we discuss some possible alternatives to mea- 
suring turbulence via the proposed magnetic fiuxes. For 
one, instead of using two SQUIDs placed at the two sides 
of the junction one could envision the use of only one 
SQUID, by changing the direction of current flow by 
merely reversing the bias. This also ensures that the 
effect of unavoidable scattering by defects/impurities on 
the magnetic fiuxes is accounted for identically for both 
possible directions of overall current flow. Another alter- 
native is to use a movable SQUID or scanning Hall- 



probe microscopy 17|. In particular, if an X-Y stage 



were added to the apparatus, so that the substrate (or 
the SQUID) were movable, one could generate images of 
the magnetic flux as a function of position. For instance, 
this scanning SQUID microscopy has been previously em- 
ployed to study the spatial configuration of vortices in 
Type II superconductors [3]. Imaging electron flow in 
this way would provide information about the electric 
current density throughout the device with the added 
benefit that the measurement would be noninvasive [2^ . 
However, we expect scanning SQUID microscopy to have 
a lower sensitivity due to the increased distance between 
the SQUID and the sample. 

One can also tune the critical current value at which 
the transition between laminar and turbulent regimes oc- 
curs by using materials with different effective masses. 
For example, the heavy-hole effective mass in p-doped 
GaAs is about 0.45me, which implies that the transi- 
tion to turbulent fiow in p-doped GaAs should occur at 

= (0.067/0.45) * /c = 0.047 AtA for the same doping 
density as the n-doped case. 
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Note that if turbulence is observed in a 2D system at a 
given current, turbulence is even more favored in 3D at 
the same current. 

Incompressibility neglects the formation of surface charges 
near the edge of the gap between the contacts 7, 19, 20|. 
These charges form dynamically during the initial transient 
of the current, and during that time create a displacement 
current. The latter would affect the initial-time magnetic 
field. However, this surface charge distribution is station- 
ary after the transient, and therefore it does not influence 
the long-time behavior of the magnetic field which is of 
interest here. 

An estimate of the flux decrease with distance can be made 
considering the magnetic field flux - through a square 
SQUID of side I - created by a current in a wire located at 
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a distance h beneath a SQUID side. It can be shown that 
the magnetic field flux is proportional to Log(l + h? /l^). 
Taking I = 200 nm, hi — 50 nm and /12 = 100 nm we 
obtain the flux ratio of 1.76. 

The flux ratio is expected to approach unity in the limit 
of zero current and infinite systems. 

Extracting the current density from the magnetic field re- 



quires solving the "inverse problem." This is nontrivial, 
but the problem has been heavily studied; for example, 
the practitoners of Magnetoencephalography (MEG) [211 ] 
have treated the inverse problem in the context of medicine 
in order to image electrical activity in the brain. 



